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Asymptotics with a positive cosmological constant: 

II. Linear fields on de Sitter space-time 
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Linearized gravitational waves in de Sitter space-time are analyzed in detail to 
obtain guidance for constructing the theory of gravitational radiation in presence of 
a positive cosmological constant in full, nonlinear general relativity. Specifically: i) 
In the exact theory, the intrinsic geometry of I is often assumed to be conformally 
flat in order to reduce the asymptotic symmetry group from Diff(X) to the de Sitter 
group. Our results show explicitly that this condition is physically unreasonable; ii) 
We obtain expressions of energy-momentum and angular momentum fluxes carried 
by gravitational waves in terms of fields defined at iii) We argue that, although 
energy of linearized gravitational waves can be arbitrarily negative in general, grav¬ 
itational waves emitted by physically reasonable sources carry positive energy; and, 
finally iv) We demonstrate that the flux formulas reduce to the familiar ones in 
Minkowski space-time in spite of the fact that the limit A —)• 0 is discontinuous 
(since, in particular, X changes its space-like character to null in the limit). 

PACS numbers: 04.70.Bw, 04.25.dg, 04.20.Cv 


I. INTRODUCTION 

A rich theory of isolated gravitating systems, developed systematically since the 1960s 
umi, lies at the foundation of a large fraction of research in general relativity with zero 
cosmological constant. Examples include the gravitational radiation theory, classical and 
quantum aspects of black holes, and several major initiatives in geometrical analysis (see, 
e.g., [5] for a summary). But observations strongly indicate that the cosmological constant 
is positive in the universe we inhabit [6] . Therefore it is important to extend the conceptual 
framework from the A = 0 case to the A > 0 regime. 

In the first paper in this series [Hj we began an exploration of this problem. Our findings 
related to gravitational waves can be summarized as follows. If one considers space-times 
which are asymptotically de Sitter in the sense introduced by Penrose |1] (more precisely, 
which satisfy Dehnition 2 in [5]) then the asymptotic symmetry group is simply Diff(X). 
Thus, with these boundary conditions, one cannot single out translations or rotations even 
asymptotically} Consequently, one cannot introduce 2-sphere charges analogous to the 
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Bondi 4-momentum at I Hi. or calculate fluxes of energy, momentum and angular mo¬ 
mentum carried away by gravitational waves [7]. We then examined a common strategy 
to reduce Difl'(X) to the de Sitter group by strengthening the boundary conditions. The 
idea is to restrict oneself to those space-times for which the intrinsic 3-metric Qab on X is 
conformally flat. This additional restriction seems natural, because the condition is satisfied 
in the familiar examples, including the Kerr-de Sitter and Friedmann-Lemaitre space-times. 
Furthermore, the 2-sphere charges at X associated with the Kerr-de Sitter time-translation 
and rotation yield the expected mass and angular momentum [5]. 

However, we showed that the additional boundary condition is equivalent to demanding 
that the magnetic part Bab of the leading order asymptotic Weyl curvature must vanish at 
X. Now, in the case of Maxwell fields on asymptotically de Sitter space-times, the analo¬ 
gous requirement would be that the magnetic held Ba should vanish at X. This requirement 
would remove half the space of solutions by fiatl By analogy, in the gravitational case, the 
strengthening of the boundary conditions appears to be physically unjustifiable. Further¬ 
more, irrespective of whether one strengthens the boundary conditions in this manner or 
not, one does not have expressions of fluxes of energy-momentum and angular momentum 
carried away by gravitational waves. Indeed, for A > 0, no gauge invariant characterization 
of gravitational waves is available in full general relativity! Thus, we have an apparent im¬ 
passe: On the one hand the Bab = 0 condition is too strong but, on the other hand, the 
boundary conditions are too weak without it (both for the gravitational radiation theory 
and quantum considerations). 

A new framework is being constructed to overcome this difficulty and address related 
issues discussed in [5]. In this paper we will complete the first step of that program by 
analyzing source-free, linearized gravitational waves in de Sitter space-time. In the A = 0 
case, the analogous analysis of linearized fields was necessary for the derivation of energy 
loss due to a time changing quadrupole moment in the weak field approximation (see, e.g., 
i)- More generally, it provided considerable intuition and important checks in the final 
construction of the theory of gravitational waves in exact general relativity [IHai9]. In 
subsequent papers we will see that same is true in asymptotically de Sitter space-times. 

The main ideas of this paper can be summarized as follows. 

We will restrict ourselves to the (future) Poincare patch of de Sitter space-time because, 
as we will see in [lO], this provides the setting that is appropriate for describing isolated 
systems in full general relativity. The subgroup of the 10-dimensional de Sitter group that 
leaves this patch invariant is 7-dimensional, consisting of 4 (de Sitter) translations and 3 ro¬ 
tations; symmetries that enable one to define the total 4-momentum and angular momentum 
carried by test fields, including linearized gravitational waves. Because of the high degree of 
symmetry of the Poincare patch, as is well-known, one can solve linearized Einstein’s equa¬ 
tion explicitly. By examining the behavior of solutions at X"*" we will explicitly show that 
the (linearized analog of the) condition Bab = 0 at X'*' removes, by hand, half the number of 
degrees of freedom associated with gravitational waves. This will confirm the expectation 
from Maxwell’s theory. 

For test matter, such as scalar. Maxwell or Yang-Mills fields, conserved quantities can be 
readily constructed using the stress-energy tensor. For linearized gravitational fields, on the 
other hand, we do not have a gauge invariant, local stress-energy tensor because in general 


states which are necessary, e.g., to systematically discuss whether the quantum evaporation of black holes 
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relativity gravity is absorbed into space-time geometry. Therefore a new strategy is needed. 
A convenient route is provided by the covariant Hamiltonian framework where the phase 
space Tcov consists of solutions to linearized Einstein’s equation. Diffeomorphisms generated 
by isometries have a well-defined action on Tcov which preserves the natural symplectic 
structure 0 on Tcov The Hamiltonians generating these canonical transformations provide 
us with formulas of energy-momentum and angular momentum carried by gravitational 
waves. We will express these quantities in terms of helds that are well-dehned on X+. 
These expressions will be needed in the derivation of the energy loss due to a time-changing 
quadrupole moment in the A > 0 case, derived in m- 

Finally, we discuss the A —>■ 0 limit. Physically one expects that in this limit energy- 
momentum and angular momentum expressions should reduce to the well-known ones for 
linear gravitational waves in Minkowski space. However, the limit is delicate because of 
conceptually important discontinuities. In particular, while is space-like for every A > 0, 
it is null for A = 0. Similarly, while the generator of every de Sitter ‘time translation’ (used 
to define de Sitter energy) is space-like in a neighborhood of X"*" for any A > 0, it is time-like 
in a neighborhood of X’*' for A = 0. Consequently, while the flux of energy at de Sitter X+ 
can be arbitrarily negative no matter how small A is, it is strictly positive in the A = 0 
case. We provide a detailed, systematic procedure to take the limit and show that the de 
Sitter fluxes do go over to the Minkowski fluxes in the limit. This procedure will be useful 
in reliably estimating the errors one makes by working in the asymptotically flat context 
rather than asymptotically de Sitter. 

The paper is organized as follows. In section we collect results on the geometry of the 
Poincare patch that will be used throughout our discussion and show that, for test Maxwell 
helds, the familiar huxes of (the de Sitter) energy-momentum and angular momentum ob¬ 
tained using the stress-energy tensor can be derived using Hamiltonian methods that do 
not refer to the stress-energy tensor. In section |IH| we study the asymptotic behavior of the 


explicit solutions to the linearized Einstein’s equation in the Poincare patch and analyze the 


consequences of the Bab = 0 condition. In section [IV] we introduce the covariant phase space 
Tcov of linearized gravitational helds in the Poincare patch, derive expressions of Hamiltoni¬ 
ans associated with the seven isometries, and express their limits to using helds that have 
well-dehned limits there. In section jVj we derive two properties of these huxes. First, we 
show that when the subtleties associated with the A —)■ 0 limit are taken into account, our 


hux expressions of section do reduce to the standard hux formulas associated with linear 
gravitational waves in Minkowski space-time. Second, we show that although gravitational 
waves in de Sitter space-time can carry arbitrarily large negative energy, for the class of 
solutions that are of direct physical interest in the investigation of isolated systems, they 
carry positive energy. 

Our conventions are as follows. Throughout we assume that the underlying space-time is 
4-dimensional and the space-time metric has signature The curvature tensors are 


dehned via: 


2VfaVq/Cc Rabc kdi 


Rac Rabc aud R Rabf] 


ab 


II. PRELIMINARIES 


This section is divided into two parts: i) symmetries of the Poincare patch; and, ii) the 
covariant phase space and conserved quantities associated with these symmetries. 
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FIG. 1: Left Panel: The Penrose diagram of a spherical isolated star in general relativity with 
A > 0. The solid diagonal line denotes the future event horizon of i~. The star and the 

radiation it emits are invisible to all observers whose world-lines are confined to the lower portion 
of the de Sitter space-time below Therefore in the discussion of this isolated system, it is 

natural to restrict oneself to the upper half. The dashed diagonal line is E~{i'^), the past event 
horizon of T*". Right Panel: The Poincare patch of de Sitter space-time of interest is the upper 
triangle, to the future of the event horizon E'^{i~) where rj = —oo. The rj = const lines denote the 
cosmological slices, i.e., flat Cauchy surfaces. 


A. The Poincare patch 

In the A = 0 case, to study isolated systems in the weak held limit, one investigates 
linearized gravitational helds in Minkowski space-time. For the A > 0 case, it may seem 
natural to replace Minkowski space with de Sitter space-time. However, because of the 
differences in causal structures of these two space-times, an important difference arises. 
Consider an isolated system -such as a single star or a binary- that is conhned to a spatially 
bounded world-tube for all times (see the left panel in Fig. [^. In this case the matter 
world-tube has future and past end-points in both A = 0 and A > 0 cases, denoted by 
However, whereas in the A = 0 case the future of i~ is the entire Minkowski space-time, 
if A > 0, it is only the future Poincare patch of de Sitter. No observer whose world-line 
is conhned to the past Poincare patch can see the isolated system or detect the radiation 
it emits. Therefore, to study this system, it suffices to restrict oneself just to the future 
Poincare patch rather than the full de Sitter space-time. Indeed, while it is difficult to 
impose the physically appropriate ‘no incoming radiation’ boundary condition at lEl. 
as we will see in [10], this condition can be naturally imposed at the cosmological horizon 
that constitutes the past boundary of this Poincare patch. Because our primary 
purpose is to develop intuition for the full, nonlinear theory, we will restrict ourselves to 
this future Poincare patch, although all our results can be readily extended to the full de 
Sitter space-time. Next, for easy comparison with the rich literature on gravitational waves 
in cosmology, we will use coordinates t], x, y, z] with x, y, z assuming their full range on 
and the conformal time y G (—cxo,0) (see the right panel in Fig. [^. Then the de Sitter 
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metric can be expressed as: 

Qabdx'^dx’’ = {1/Hr]f {- dr]^ + dx"^) =: {1/Hr]f gabdx^dx\ ( 2 . 1 ) 

where H := =: 1/^ is the Hubble parameter, the inverse of the cosmological radius 

^ While these coordinates are extremely convenient in the detailed calculations of 
gravitational perturbations, it is obvious that they are ill-suited for taking the limit A —)■ 0. 
To take this limit, it is simplest to use proper time f, which is related to the conformal time 
r] via Hrj = —In terms of t, the de Sitter metric becomes 

Qabdx^dx^ = -df + 2 ( 2 . 2 ) 


and it is manifest that the metric coefficients go to those of Minkowski metric coefficients as 
A goes to zero. Therefore, to compare geometric structures in de Sitter space-time to those 
in Minkowski space, it is important to use the differential structure induced on the Poincare 
patch by {t,x), and not by {r],x)\ 

Locally, of course, this metric admits 10 (de Sitter) Killing helds. However, since the 
Poincare patch is only a part of the de Sitter space-time, only those isometries are permissible 
that map this patch to itself. Therefore, we now have to restrict ourselves only to those 
Killing fields that are tangential to its boundary in the full de Sitter space-time. 

As discussed in detail in section 4.C.2 of |5], these Killing fields constitute a 7-dimensional 
family. We have 3 spatial translations and 3 spatial rotations tangential to each 
T] = const slice, generating the Euclidean group. In addition, there is a 7th Killing field 


T 


-H 


d d d d 

^ ^w 

or] ox oy oz 


(2.3) 


We will refer to T“ as time translation because: i) it is the limit of the time translation 
Killing field in the Schwarzschild-de Sitter space-time as the mass goes to zero, and, ii) in 
the {t, x) coordinates, it reduces to a time-translation in Minkowski space-time as A —)• 0.^ 
The commutation relations between these seven Killing helds are given by: 


[^5 ^h)] = ^ ’ [^5 -^h)] = O’ [^(b’ [-^(b’ -^0')] = ■ 

(2.4) 

(Note that the time translation does not commute with space-translations.) We will denote 
this 7-dimensional Lie-algebra of symmetries of the Poincare patch by gpoin and the Lie 
group it generates by C5poin-^ Finally, in the standard conformal completion of the Poincare 


^ These coordinates -as well as the coordinates {t,x) discussed below- have the disadvantage that they 
do not cover the past boundary of our Poincare patch, i.e., the event horizon of i~. But this 

limitation will not affect our considerations. 

^ Limit A —)► 0 of illustrates the importance of using the correct differential structure to take this limit. 
Had we used the differential structure provided by {r],x) we would have concluded from (2.3) that T® 
vanishes in the limit. But this procedure would have been incorrect because the metric (jab diverges in this 
limit (although it reduces to the well-defined Minkowski metric if the limit is taken using the differential 
structure induced by {t,x)). Note, incidentally, that T“ is sometimes referred to as ‘dilation’ because it 
is the conformal-Killing vector field representing a dilation with respect to the flat metric gab- 
^ This is the group that leaves the point i~ on I~ of de Sitter space-time invariant. As A —>■ 0, ©poin 
reduces to a well defined seven dimensional subgroup of the Poincare group; the limit carries the memory 
of the preferred t = const slicing. 
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patch, X"*" has topology and this 7-dimensional gronp preserves the completeness of the 
allowed class of metrics on X+ [5] . 


B. Maxwell fields in de Sitter space-time 

As is well-known, each Killing symmetry 77“ leads to a conserved qnantity. For matter 
fields -such as the Maxwell field Fat- the standard procedure is to use the stress-energy 
tensor Tab = ^{FamFbng"^'' - {l/F)gabFcdFmng‘''^g‘^'^) ■ The conserved quantity associated 
with a Killing field 77“ is given by 


Fk= [ r,b77“n'dVs 

Jt. 


( 2 . 5 ) 


where the integral is taken over any Cauchy surface S with unit normal n“. Fk may be 
regarded as the ‘flux’ of the conserved quantity across S. 

However, for the linearized gravitational field, we do not have a gauge invariant, locally 
defined stress-energy tensor. We will now show that, in the Maxwell theory, the expression 
(2.5) of Fk can also be obtained using a covariant phase space framework without having to 
refer to the stress-energy tensor. In section we will use this alternate method to calculate 
conserved quantities for the linearized gravitational field. 

Consider a globally hyperbolic space-time, (M,gab) with a Killing field 77“. Denote 
by 7^“^^ the space of all suitably regular, source-free solutions Fab to Maxwell equations 
V[aXfec] = 0 and cFab = 0. Starting from the Maxwell Lagrangian, one can show that 
naturally endowed with a symplectic structure (i.e., a closed, non-degenerate 2-form) 

^Max- 

WMax(F, = ^ f [FabAa - FabA^ ] u' dVs . (2.6) 

i/ s 


Here F and X are any two solutions to Maxwell equations, Aa is any vector potential for 
Fab (he.. Fab = 2 V[aA;,]) and S is again any Cauchy surface. Using Maxwell equations (and 
the suitable fall-off implicit in the regularity condition) it is easy to verify that the right 
side is independent of the choice of the Cauchy surface S and is gauge invariant. The pair 
, i^Max) is the Maxwell covariant phase space. Each Killing field 77“ on Ai naturally 
defines a vector field /C on via: IC\f= 6kF := CKFab. Not surprisingly, the flow 

generated by K. on F^^ preserves the symplectic structure WMax, i-e., defines a 1-parameter 
family of canonical transformations on (F^y’^, (UMax)- The Hamiltonian generating this flow 
is a function TLk on given by: 


F-K '■= —2'^Max(X, CkF) . 


(2.7) 


For any Killing field 77“ one can verify that Fk defined in (2.7) equals Fk defined in (2.5) 


(For details on the covariant phase space of fields, including general relativity, see, e.g., [I3].) 

Let US illustrate this result for the Killing fields in the Poincare patch. Let ns first set 
77“ = 5*“, where 5*“ stands for any one of the 6 Killing fields and 7?“^^, tangential to the 
space-like slices S given by = const. Then, we have 


Fs=]- [ (FamFbn^^S'^n’^) dVs = ^ 

dvr Js 47r 


{eabcE'^B'^S’^) dVs 


( 2 . 8 ) 
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where Ea := FabU^ and Ba '■= *Fabn^ are the electric and magnetic parts of the Maxwell 
held, and eabc the alternating tensor on the slice S. Thns, as one would expect, Fs is the 
hux of the S'-component of the Poynting vector eabcE’^B^ across S. Next, let us consider the 


Hamiltonian (2.7) generated by S\ 


^5 = -^ / r'^[{CsFab)A, - 

= ^ f^F^EsAb) n“ dVs = ^ {eabcE^B^S’^) dVs (2.9) 

where in the hrst step we have integrated by parts and in the second step used Cartan 
identity and the Maxwell equation DaE°^ = 0. Thus, using the covariant phase space we can 
recover the conserved quantity Fs as the Hamiltonian Fs dehned by the Killing symmetry 
Because of conformal invariance of Maxwell equations, we can easily take the limit as E 
approaches X’*' and express the conserved hux as an integral over . The expression (2.8) 


brings out the fact that if the magnetic held vanishes at X+, then that electromagnetic wave 
carries no angular momentum or linear momentum. 

For the time translation T“, the argument establishing the equality of Fk and XCk is the 
same but the calculation is a little more involved because T“ has components both along 


and orthogonal to the cosmological slices (see Eq.(2.3)). We hnd: 


.Ft = Wt = ^ / [(BaBi + r" + 2e^tcE'‘B‘r] dV t. 


( 2 . 10 ) 


In the limit as E approaches X^, T“ becomes tangential to X^ (since r/ = 0 at X^) and 
vanishes. Therefore the expression of the conserved energy reduces to an integral of the 
component of the Poynting vector along T“: 


Ft = 'Ht = ^ [ {eabcE^B^T^) dVx+ (2.11) 

dvr Jx+ 

where the electric and magnetic helds and the alternating tensor are calculated using any 
conformally rescaled metric that is regular at X+ (e.g., gab)- This expression brings out two 
interesting facts. First, in de Sitter space-time while the energy carried by electromagnetic 
waves is conserved as in Minkowski space-time, now it can be negative and is unbounded 
below. Second, if we restrict ourselves to Maxwell helds whose magnetic held vanishes at 
X+, then those electromagnetic helds carry no energy either. Note that the second result 
is specihc to X"*": If the magnetic held vanishes on a cosmological slice rj = const ^ 0, 
the energy of that Maxwell held does not vanish unless the Maxwell held itself vanishes 
identically. The 3-momentum and the angular momentum, on the other hand do vanish. 

To summarize, for Maxwell helds, the conserved quantities associated with Killing helds 
in the Poincare patch can be recovered as Hamiltonians on the covariant phase space, 
without any reference to the stress-energy tensor. Also, because all Killing helds K°' on 
de Sitter space-time are tangential to X~^ -and hence space-like- one can express every 
conserved quantity Fk as an integral across X+ of the component of the Poynting vector 
along 77“. This expression brings out the fact that if we were to require that the magnetic 
held vanish at X+, we would be left with electromagnetic waves that carry no 3-momentum 
or angular momentum, nor energy defined by de Sitter isometries! 





III. LINEARIZED GRAVITATIONAL FIELDS 


As in section II A, we will use the (r/, x) chart and the form (2.1) of the de Sitter metric 


gab in the Poincare patch. The perturbed metric will be denoted by gab, 


9ab 9ab T ^ 'lab 


(3.1) 


where e is the smallness parameter and 'jab denotes the hrst order perturbation. Then, in 
the Lorentz and radiation gauge, i.e., when the gauge freedom is exhausted by requiring 
that 7 afe satisfy 

= 0; 7 ^ 67 “ = 0; and 7 ^ 6 ^“'’ = 0, (3.2) 

the linearized Einstein’s equation simplihes to 

n^ab-2H^^ab = 0. (3.3) 


(Here is a vector field normal to the cosmological slices with r]^da = d/dr].) Following a 
common strategy in the cosmology literature, it is convenient to rewrite (3.1) as 


9ab — ^ (^) i,9ab H” ^ 


{^9ab “ 1 “ ^ f^ab) 


(3.4) 


(HvY 

since calculations are simpler in terms of the mathematical held hab than in terms of the 
physical perturbation 7 ^ 6 . Indeed, the gauge conditions (3.2) can now be written using the 
background flat geometry of gab- 


^""hab = 0; habT]'' = 0; and habg'"’' = 0 , 
and the linearized Einstein’s equation becomes 


Uhab-2—h'ab = □ hafe + - = 0 , 

a rj 


(3.6) 


(3.6) 


where Jiab- Note that the gauge conditions and linearized Einstein’s equation sat- 

ished by hab are the same as those satished by the linearized gravitational helds in Minkowski 
space-time in absence of a cosmological constant except for the extra term ( 2 / 7 ) in the 
linearized Einstein’s equation. In the {t, x) differentiable structure that is well-suited to take 
the limit A —)■ 0, the extra term {fl/tfldr^hab = —2Hdthab goes to zero, just as one would 
expect. 

As in the case of linearized helds in Minkowski space-time, it is simplest to hnd explicit 
solutions using a Fourier transform: 


Ki{x,r]) = j e“' 


(3.7) 


where (s) labels the two helicity states and are the polarization tensors, satisfying 


( 7=1 
(«)/ 


i1|(^=0i cl7fc)fc“ = 0; cl7fc)r = 0; 




(s)/ 


(ci?(*;))* = 7 ((-*;); 


— 


^ab 




ac 9,bd 


(3.8) 
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Here, and in what follows, qab is the fixed spatial Enclidean metric on the cosmological 
slices, tailored to the co-moving coordinates x, and * denotes complex conjngation. The two 
fnnctions h^^\ri) captnre the gange invariant information -the transverse traceless modes- 
of the linearized gravitational held. Since hab{x, rj) are real helds, it follows that 


(/!<■>(,))* = . 


( 3 . 9 ) 


The held eqnation (3.6) implies that the satisfy the ordinary diherential eqnation (ODE); 


(4“')" - -(ft'-i)' + ehf-' = 0, 


r] 


(3.10) 


where the prime denotes diherentiation with respect to rj, and k'^ = k ■ k. The second order 
ODE ( |3.10[ ) can be readily solved to obtain the general solntion 

h^^\v) = (~2if) [r] cos{kri) — {1 /k) sin{kri)) — [j]sin{kri) + {1 /k) cos(fc77))] (3.11) 


where and are arbitrary coefficients (in the Schwartz space), determined by the 
initial data of the solntion. (These coefficients can also depend on A. We did not make this 
dependence explicit becanse in the main text we work with a hxed valne of A.) Snbstitnting 
(3.11) in (3.7) we obtain the general solntions hab representing hrst order pertnrbations. 

Next, let ns discnss cnrvatnre. Since the Weyl tensor of de Sitter space-time vanishes, 
the hrst order pertnrbations ^^'^Eab and ^^^Bab of the electric and magnetic parts of the Weyl 
cnrvatnre are gange invariant and can be expressed directly in terms of the solntions hj^{ri) 
in (3.11). To hnd these expressions, we hrst note that, in exact general relativity, the electric 
and magnetic parts are related to the hrst and second fnndamental forms Qab and Kab on 
any space-like snrface via 


Eab = TZab - K^K^b + KKab - + qabq^"){^Rmn - (1/6) ^Rgmn) 

Bab = e(^a^^ D\m\K\n\b) (3.12) 


where D, eabc and TZab are the derivative operator, alternating tensor and the Ricci cnrvatnre 
of the 3-metric qab, and '^Rab is the Ricci cnrvatnre of the space-time metric gab- 

It is straightforward to linearize these eqnations nsing the cosmological foliation on the 
de Sitter backgronnd. Calcnlations are simplihed by noting that: (i) Eab and Bab are con¬ 
formally invariant, and, (ii) a convenient conformal completion of de Sitter is provided by 
choosing the conformal factor D = — Hr], so that the conformal metric kEcjab that is well 
behaved at is just the Minkowski metric gab in the {rj, x) chart. Therefore, in ehect, 
linearization can be carried ont using this hat background metric. The perturbed electric 
and magnetic parts of the Weyl tensor can be expressed using hab and geometric structures 
associated with the hat 3-metric qab on each cosmological slice: 

(^)Eab = -l{D^hab + -h'a,), aud ^^^Bab = ha^^Dim\h\^\by (3.13) 

Zi Tj Z 

Recall that the boundary conditions at imply that the Weyl curvature of an asymptoti¬ 
cally de Sitter metric must vanish at [HE]. Therefore, the hrst order perturbations ^^"^Eab 
and ^^^Bab of Weyl curvature also vanish at X+ and 

:= fi-‘ and := 


( 3 . 14 ) 
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admit smooth limits there. We will refer to £ab and the Bab as the perturbed eleetrie and 
magnetie parts of the Weyl eurvature as a short hand since it is these quantities that will 
feature in most of our discussion. Using explicit solutions (3.11) it is easy to verify that they 
do indeed admit smooth limits to X+: 


Babij^i h) |f?=0 


Bab{x, 7]) 1^=0 = 


(27r)3 

d^k 

(27r)3 




(s)=l 


^ erf(k)e 


ik-x 


(3.15) 


(s)=l 


where (kn/k) is the ‘dual’ of the polarization tensor. These formulas bring 


(s) is) 

out the meaning of the coefficients and that feature in the expression (|3.1ip of a 


(s') 

general solution to the linearized equations. ^ directly determines the electric part of the 

perturbed Weyl tensor at and the magnetic part at X^. 

It is therefore clear that the perturbed Weyl tensor has no magnetic part Bab{x, rf) at X+ 
if and only if the solution hab has the form 


h^f\ri) = {—2H) (t] cos kr] — (l/k) s'mkr])] 


(3.16) 


everywhere, obtained by setting = 0 in (|3.11|). Thus, the condition that the magnetic 


part vanish at X~*" -or, that conformal flatness of the 3-metric at X+ be preserved to 
hrst order- removes, by fiat, half the degrees of freedom from consideration. The hrst 


expectation based on Maxwell fields is explicitly borne out. In section IV we will show 
that the second expectation is also borne out: the remaining gravitational waves, that do 
preserve conformal flatness to hrst order, carry no energy, momentum or angular momentum. 


i?emarks: 


(i) The explicit solution (3.11) shows that, as one approaches X'^ (i.e. as p ^ 0), the term 
associated with Ui ^ vanishes while the term associated with Si survives. In the cosmology 
literature, the hrst is referred to as the ‘decaying mode’ and the second as the ‘growing mode’. 
Thus, the requirement that the magnetic part of the perturbed Weyl curvature vanish at 
removes by hat the growing mode and leaves only the decaying mode. These perturbations 
hab vanish at . 


(ii) Let us return to the linearized Einstein’s equation (3.6) satished by hab- While one can 
think of hab as a held propagating on the Minkowski metric gab, because of the additional 
term {2/ri)h'^f„ the propagation is not sharp; there is a ‘tail term’. On the other hand, 
the linearized Weyl tensor satishes conformally invariant equations. Its propagation does 
not have a tail term. Interestingly, the same is true of the time derivative of the metric 
perturbation: One can verify that it satishes the conformally invariant equation, (□ — 
{fR/d))h'^^ = 0. Equivalently, since gab = i^/9ab, h follows that O [{1 /rj) h'^i^] = 0. 
Therefore it follows that the propagation of {l/ri)h'^f, on (M,gab), and hence of h'^f^ on 
{A4, (Jab), is in fact sharp, without any tail terms. This fact has an interesting implication 
in the discussion of the quadrupole formula m- 
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IV. THE HAMILTONIAN FRAMEWORK 

This section is divided into two parts. In the hrst, we constrnct the covariant phase 
space of sonrce-free, linearized gravitational helds on the de Sitter backgronnd. In the 
second, we obtain expressions of energy, momentnm and angnlar momentum carried by 
gravitational waves by computing the Hamiltonians corresponding to the seven Killing helds 
on the Poincare patch. 


A. The covariant phase space 


For linearized gravitational helds, the covariant phase space Fcov can be taken to be the 


space of solutions jat to the equations (3.2) and (3.3). For simplicity, we will assume that the 


solutions of interest have initial data in the Schwartz space of rapidly decreasing, smooth 
helds, although these conditions can be weakened considerably. The standard procedure 
e.g. [13]) endows Fcov with a symplectic structure u. Restricted to the cosmological 


see, 


slices S (given hy rj = const), it becomes: 


uj{y,y) = uj{h, h) := 


a\v) 


d^x{habllcd - Kbhcd) 


(4.1) 


where hab is related to the physical metric perturbation jab via jab = o? hab (see Eq. (3.4)) 
and K = SttG. It is easy to verify that (3.5) and (3.6) imply that the integral is independent 


of the T] = const slice on which it is evaluated. This form of the symplectic structure is useful 


in calculations within the Poincare patch. Furthermore, as we will see in section IV B, it is 
well-adapted for taking the limit A —)■ 0. 

In the cosmology literature, one often works with the functions h^^\ri) dehned in (3.7) 
and their Fourier transforms 


(j)^^\x,ri) : = 




(^ k 




iv) e 


ik-x 


( 4 . 2 ) 


in place of the tensor helds yab or hab- (The factor of \/4fi; is introduced to endow with 
the standard dimensions of a scalar held, so that the scalar and tensor perturbations can be 
treated in a completely parallel manner. See, e.g., section 3.D of [Il|.) These are referred 
to as the two tensor modes. It is straightforward to verify that these helds satisfy the wave 
equation in de Sitter space-time 

n(f)^^\x,r]) = 0 . (4.3) 

Thus, each tensor mode 0^^^ of the linearized gravitational held satishes just the massless 
Klein-Gordon equation. It is clear that, given hxed polarization tensors t^ab(^)i there is a 
natural isomorphism between the functions 0^'^^ {x, r]) and solutions hab{x, rj) to the linearized 
Einstein equation (3.6) and gauge conditions (3.5). It is easy to check that the symplectic 


structure (4.1) on Fcov translates to the standard symplectic structure on the covariant phase 
space F^^ consisting of pairs of solutions 0 = {(p^^\x,r])} to the Klein-Gordon equation: 




' P)=i 


(4.4) 
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This form of the symplectic structure is useful to compute expressions of fluxes of energy- 
momentum and angular momentum that are adapted to the ‘tensor modes’ used in the 
cosmological perturbation theory. 


However, expressions (4.1) and (4.4) of the symplectic structure have one drawback: 
because of the multiplicative factor aF(j]) = {1 /H'^t]'^), they are not well-suited to take the 
limit to (where t] = 0). While, the limit itself is well defined because the symplectic 
structure is independent of rj, to express physical results -e.g. the formula of energy- in 
terms of helds that are well defined at X'*', one has to be extremely careful in keeping track 
of terms in the integrand which tend to zero at the appropriate rate to compensate for 
the apparent blow up as l/r/^ due to the pre-factor in front of the integral. Also, these 


expressions are not gauge invariant as they use specihc gauge conditions (3.5). To overcome 


these limitations, it is convenient to recast the expression (4.1) using the relation between 


the perturbed electric part of the Weyl tensor Sab and the metric perturbation. 


t/) = — (h'^+riD^h^i) , 


(4.5) 


that holds on any cosmological slice. Substituting for h'ah in terms of Sab and simplifying by 
performing integrations by parts, we obtain: 


u;{h, h) = 


2Hk 




( 4 . 6 ) 


We will use both expressions, (4.1) and (4.6), of the symplectic structure on Tcov in our 


discus sion of the conserved fluxes associated with the 7 Killing vectors. (The equivalent 
form (4.4) in terms of the Klein-Gordon helds turns out not to be as useful in providing 


hints for the full, nonlinear theory.) 

We will conclude this discussion by pointing out several consequences that follow im¬ 


mediately from the form (4.6) of the symplectic structure. First, it is transparent that 
{1/2Hk)S'^^ can be regarded as the momentum that is canonically conjugate to the metric 


perturbation hab- Second, as we saw in section |III[ the perturbations hab as well as the 
perturbed electric part of the Weyl tensor Sab admit well dehned limits to X+. Therefore, 


one can take the limit S —simply by evaluating the integral (4.6) on X+. This feature 


will facilitate our task of expressing energy, momentum and angular momentum in terms 
of asymptotic helds at X+. In turn, these expressions will be directly useful in [TT] to ob¬ 
tain a formula for the energy emitted by a time changing quadrupole, and establishing its 
positivity. The third and more important feature is gauge invariance. Note hrst that Sab 
by itself is gauge invariant, it is tangential to the cosmological slices, and it is divergence- 


and trace-free. This fact enables us to drop the gauge hxing conditions (3.5) and consider 


general perturbations. For, if either 'fab (or, 'fab) is a pure gauge held -i.e. of the form 


V(a'Cb) for a space-time vector held properties of Sab ensure that the expression (4.6) 
of u{h, h) vanishes identically. Thus, the passage from to Sab using (4.5) has provided 
us with a manifestly gauge invariant expression (4.6) of the symplectic structure. Finally, 


using the explicit solutions (3.11), we can re-express the symplectic structure in terms of 


the coefficients and 


u(h, h) = — 

K 


/ (Cl E - (A'l* Ef) . ( 4 . 7 ) 

^ > (s)=l 
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where, as before, * denotes complex-conjugation. Consequently, the pull-back of the sym- 
plectic structure to the subspace of Fcov on which Bab vanishes on 1+ -or, alternatively, 
on which Sab vanishes on X’*’- is identically zero. These subspaces are among the maximal 
Lagrangian subspaces of Tcov In this respect the situation is again completely parallel to 
that in the Maxwell theory. 


B. 3-momentum, angular momentum and energy carried by gravitational waves 


We can now calculate the Hamiltonians on Tcov corresponding to the seven Killing helds 
on the Poincare patch. Recall from (3.4) that the physical metric perturbation is 'jab = o?hab 


and it satisfies the gauge conditions (3.2) and linearized Einstein’s equation (3.3) that refer 

'-' — '-' (K^ 

only to the background de Sitter metric gab- Therefore, if '^ab G Tcov, then so is := 
S-Klab, for any Killing held iP" of Qab- From the dehnition (3.4) of hab, it follows that 




= a {Ck hab + 2(a Cxa) hab) =■ a h\^ 


(4.8) 


with a = —l/{Hri). As in the Maxwell case, the isometries generated by each of the seven 
Killing helds in 0poin provide a 1-parameter family of canonical transformations on Tcov 
From general results on the covariant phase space na it follows that the corresponding 
Hamiltonian is again given by 


Hi- := -^w(7, 7™) = ft™). 


(4.9) 


Recall from section that if Bab = 0 at X+, then hab also vanishes there. In this case, 
then, we have l-Ls = 0. Thus, although there do exist linearized gravitational waves that 
retain conformal hatness of the induced geometry at to hrst order, they carry no energy, 
3-momentum or angular momentum. 


We will now compute the Hamiltonians (4.9) for the seven Killing helds in gp, 


1. 3-momentum and angular momentum 


As in the case of Maxwell helds discussed in section H B 
for the 3 spatial translations and the 3 rotations 


the calculations are identical 


Let us therefore again denote by 
S°' any of these six Killing helds and calculate the 3-momentum or angular momentum XLs, 
and then discuss energy XLt separately. For these six Killing helds, we have h\f = Cs hab 


since these helds are all tangential to the p = const surfaces. Furthermore, from (4.5) it 


(S) 

follows that the corresponding perturbed electric part of the Weyl tensor, is given by: 


p('S’) _ 

^ab 


2Hri‘^ 


[Crf{Cshab) +3 D'^iCshab)) = Cs Sab 


Therefore (4.9) becomes: 


(4.10) 


Us = -^uj{h, h(^)) = 






s 

ac °cd 


2Hk 


(f„j CsKJ) q“‘q' 


(4.11) 
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where, in the second step we have integrated by parts. Thus, the expressions of 3-nionientuni 


and angular momentum mirror those in the Maxwell theory. Since the integrand in (4.11) 


refers only to the helds hab, £-ab and the metric all of which have smooth limits to X’*', 


to take the limit E —)■ X"*" we just have to evaluate (4.11) on X+. 


Finally, let us consider the limit A —)■ 0 of "Hg. Since the Hubble parameter H tends to 


zero in this limit, from the form of (4.11), the limit seems divergent at hrst sight. However, 


this conclusion is incorrect because helds in the integrand also depend on H. Let us therefore 
analyze the limit more carefully. As explained in section [TT| to take this limit, we should use 
the differential structure induced by the chart (f, x) on the Poincare patch (and not by the 
chart (? 7 , x)). Then, for the background geometry, we hnd that as A —)■ 0, we have 


— 


Hr] = —e 


gab -t Vab = -dot dbt + daX dbX + daV BbU + BaZ BbZ] 

d d 


— 1 ; 


— —^ — = Hf) ' 
dr] dt 




(4.12) 


Note that the Minkowski metric r]ab in (4.12) is distinct from the Minkowski metric gab in 
(3.4). In the Poincare patch, each of the Cartesian coordinates, {t, x) of r]ab takes the full 


range of values, (—oo, cx)) (whereas r], of gab only runs from (— cxd, 0)). A second important 
point for the limit is that the {t, x) chart does not cover X+ (where t = oo). Therefore, to 
take the A —)■ 0 limit, we are led to evaluate the symplectic structure and Hamiltonians Hs 
on a cosmological slice corresponding to a hnite constant value of t. 



show that hab satishes the Lorentz and radiation gauge, as well as the linearized Einstein’s 
equation with respect to the Minkowski metric r]ab'. 


d'^hab = 0; habV = 0; habf = 0; and d’^d^ hab = 0 . 


(4.13) 


Denote the space of solutions hab to these equations by fcov It is straightforward to verify 
that in the limit A —)■ 0, the symplectic structure oj on Pcov goes over to the standard 
symplectic structure u on Pcov: 


<^Ch, h) = — d^x{habdtChad) “ dtCKb) Kd) Q'"" , 


(4.14) 


where the integral is taken on a t = const slice. Finally, the limit Hs the Hamiltonian 
Hs = (-1/2) oj{h, h(^)) is given by: 


Hs = — ^sh) 

= ^ / d?x{dt(hab) CsKd) . 


(4.15) 


This is precisely the expression of the linear and angular momentum of linearized gravita¬ 
tional waves in Minkowski space-time. Thus, although the procedure of taking the limit 


A —)■ 0 is rather subtle, the de Sitter 3-momentum and angular momentum (4.11) do reduce 
to the standard conserved fluxes in Minkowski space-time. 
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Remark: 

While taking the limit, we assumed the existence of a family hab{R), satisfying the gauge 
conditions (3.5) and linearized Einstein equation (3.6) for each A, that admits a smooth 
limit hab on the Poincare patch as A —)■ 0. An explicit example of such a family is provided 
by setting in (3.11) 


(A) = Cf sin(k/H) + Df cMk/H); 
b'‘>(A) = -rf'* cos(4/B) + sm(k/H ). 


(4.16) 


A careful calculation shows that in the limit A —)■ 0 the field h^^\A) of (3.7) goes over to 

sin kt + 2D^^\k) cos kt 

and the limit hab of the family hab{A) is given by 


( 4 . 17 ) 


hab{x,t) ^ 




{2n 


'-ab 


(k) e* 


(4.18) 


W=i 


where — i C^. Thus, a general solution to the linearized Einstein’s equation in 

Minkowski space in the transverse, traceless, radiation gauge can be obtained as a limit of 
this family hab{A). 


2. Energy 


Next, let us consider the energy RLt defined by the time translation T“ of (2.3). In 


this case, the calculation is not as straightforward because: (i) the vector field T“ is not 
tangential to the cosmological slices rj = const except at X+; (ii) has an extra term 
relative to 

= Ct hab + 277 hab ; (4.19) 

(rp\ 

and, (iii) a detailed calculation shows that also has an extra term: 

= Ct Sab-HSab. (4.20) 

Once these differences are taken into account, the conserved energy-flux 77r across S can 


be calculated using (4.9). We have: 

1 


77t = 


d^x (habCrScd — ScdSrhab — ^HhabScd) 


AHk 


(4.21) 


where, as before, S is any cosmological slice. However, since T is not tangential to S, on a 
general cosmological slice we cannot integrate by parts as we did for Ris- Again the limit to 
is straightforward since all fields in the integrand have smooth limits to X+. Furthermore, 


in the limit T“ becomes tangential to X+ enabling us to simplify (4.21) further: 

1 


RLt = 


2Hk Jx+ 


d^x Scd {Srhab + 277 hab) 


77 


/ (05 ^ E (4’^‘(A’)*+2^’ (Ab*). 


(4.22) 

(4.23) 
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where in the second step we have used the explicit solutions (3.11) for hab in terms of Fourier 


modes. Since Bab = 0 if and only if = 0, the last expression makes it explicit that if 
a gravitational wave does not change conformal flatness of the intrinsic geometry at to 


hrst order, it does not carry energy. Finally, the expression (4.9) of Bk is linear in for 


all Killing helds. Therefore, Hxt = XBr for all real numbers A. For linearized gravitational 
waves on Minkowski space-time, energy is positive dehnite and vanishes if and only if the 
perturbation is pure gauge. On the de Sitter space-time, the conserved energy Bt can have 
either sign and we have an inhnite dimensional subspace of the physical, transverse-traceless 


modes for which the energy vanishes. From (4.22) it is clear that energy also vanishes if Sab 
vanishes on X+. (The other possibility, Crhab = —2Hhab on X+ is not realized because such 
perturbations would not be in the Schwartz space on of the Poincare patch, which is 
topologically |5].) 

Finally, let us consider the limit A —)• 0 of the conserved energy Bt- For reasons given 
in section IV B 1, we have to use the differential structure induced by the coordinates (t, x) 


and work with a cosmological slice in the Poincare patch with rj ^ 0. Let us again suppose 
that we have a 1-parameter family of pert urba tions hab{A.) that satisfy the gauge conditions 
3.5) and the linearized Einstein equation (3.6), and admit a smooth limit hab as A —»■ 0. As 


discussed above, hab is a metric perturbation on the Minkowski metric rjab, satisfying (4.13) 
The limit Bt of the Hamiltonian Bt = (~l/2) uj{h, is given by: 


Bt = j d^x{habdthcd - (5tha6)( 

= ^ I d^x{{dthab){dthad) + {DMiDahad)r^ 


(4.24) 


where in the second step we have used (4.13) and integrated by parts. This is precisely 


the conserved energy flux of the linearized gravitational held hab in Minkowski space-time. 


Thus, our energy expression (4.21) for linearized gravitational helds in de Sitter space-time 


does have the expected limit as A —)■ 0. Note that the limit is quite subtle and discontinuous: 
While Bt can be negative and arbitrarily large, no matter how small the positive A is, in 
the limit A —)■ 0 we obtain Bt which is positive definite! Geometrically, this occurs because 
while the Killing held T“ of de Sitter metric Qab is space-like in the ‘upper half of the 
Poincare patch’ for every A > 0, its limit, the Killing held of rjab, is time-like everywhere. 


Remarks : 

(i) In the cosmological literature, the discussion of ‘energy’ often refers to the Hamiltonians 
Bri or Bt that generate evolution along the conformal time rj or proper time t. Since and 
are not Killing helds, these Hamiltonians are not conserved. Thus, they are unrelated to 
the conserved energy Bt discussed above and are not the analogs of the standard notion of 
energy in Minkowski space-time, used in the gravitational radiation theory. 


(ii) As discussed in section IV A, in cosmology one often encodes the metric perturbations 
7ab(xh) fhe two ‘tensor modes’ r;) that satisfy the Klein-Gordon equation with 

respect to the de Sitter metric gab- On the Klein-Gordon phase space Fcw) ^^e isometry 
generated by any Killing held again dehnes a 1-parameter family of transformations that 
preserve the symplectic structure cukg- As one would expect, the corresponding Hamiltonians 
agree with the Bk obtained above for all seven Killing helds. That is, our energy-momentum 
and angular momentum expressions Bt and Bs hold both for the metric perturbations jab 
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satisfying (3.2) and (3.3) and the ‘tensor modes’ 0*^*^ satisfying the wave equation (4.3) in 
the Poincare patch. 

(iii) Finally, we note that the explicit solutions (3.11) are widely used in the cosmological 
literature on linearized gravitational waves. However, the primary interest there is on the 
effect of these gravitational waves on the polarization of the CMB electromagnetic waves. 
To our knowledge, this literature does not contain the analysis of the asymptotic behavior 
of these perturbations at I^, or the implications of the assumption that the perturbations 
preserve conformal flatness of to linear order, nor a discussion on the isometry group 
0Poin that preserves the Poincare patch, or the associated conserved fluxes XLk given above. 


V. DISCUSSION 


In the A = 0 case, there is a well-developed theory of isolated systems and gravitational 
radiation in full, nonlinear general relativity that has played a dominant role in a number 
of areas of gravitational science. In the hrst paper [S] in this series, we showed that there 
are signihcant conceptual obstacles in extending this theory to allow a positive cosmological 
constant, however small, because the limit A —)■ 0 is discontinuous. In particular, whereas 
X is space-like, no matter how small A is, it is null when A vanishes. If A were zero and 
the accelerated expansion of the universe is caused by some matter held rather than a 
cosmological constant, that held would not have the asymptotic fall-oh we are familiar with 
in the A = 0 case, and space-time curvature far away from the sources would be similar to 
that in asymptotically de Sitter space-times. Therefore, difficulties discussed in [3] would 
persist also in the A = 0 case if the observed accelerated expansion continues to inhnite 
future. To overcome these obstacles, one needs a new framework. In this paper we completed 
the hrst step to this goal by discussing linear gravitational waves in de Sitter space-time. 

Motivated by considerations of isolated systems discussed in section IIA , we focused 
on the upper Poincare patch of de Sitter space-time. Isometries generated by 7 of the 
10 de Sitter Killing helds leave this patch invariant. This group 0poin is generated by 3 
space-translations and 3 rotations that are tangential to the cosmological slices and a time 
translation that is transversal to them. Therefore, one expects well dehned notions of linear 
and angular momentum, and energy, associated with any physical held on the Poincare 


patch. We showed in section IIB that, in the case of Maxwell helds, these ‘conserved huxes’ 
arise as the Hamiltonians generating canonical transformations induced by the action of 
Killing helds on the covariant phase space P^^. Furthermore, in the A = 0 case, the 
Hamiltonian framework has been used very ehectively also for gravitational waves in full, 
nonlinear general relativity: It leads to hux integrals corresponding to the Bondi-Metzner- 
Sachs (BMS) asymptotic symmetries |^. Therefore, it is natural to use this strategy also in 
the A > 0 case. 

Since the covariant phase space consists of solutions to the held equations, in section 


TIT] we discussed the asymptotic properties of solutions to linearized Einstein’s equation in 
de Sitter space-time. In section IV we constructed the covariant phase space (Fcov,!^) of 
these linear gravitational waves. Each of the 7 Killing helds naturally dehnes how on 
Fcov that preserves the symplectic structure uj thereon, and thus dehnes a Hamiltonian XLk- 


These Hamiltonians provided us with the expressions (4.11) and (4.21) of huxes of energy- 


momentum and angular momentum carried by gravitational waves. Furthermore, we could 
express these conserved huxes in terms of helds dehned on X'*'. 

These results have a number of interesting features. First, to make the full nonlinear 
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theory manageable in the A > 0 case, at hrst it seems natural to strengthen the boundary 
conditions by requiring that the intrinsic geometry of be conformally flat, as in de Sitter 
space-time. However, almost 30 years ago Friedrich showed that the freely specihable data 
at X“ consists, up to arbitrary conformal rescalings, of a freely specihable Riemannian 
metric and a trace-free, symmetric tensor held of valence two, which satishes a divergence 
equation [19]. Therefore, by applying those results to X+ (in place of X“), it follows that 
demanding conformal hatness of the metric at X’*' removes by hand part of this free data. 
In the linearized approximation, we could sharpen the implication of this condition. First, 
because the perturbed electric and magnetic parts of the Weyl curvature are gauge invariant, 
we can discuss the physical or true degrees of freedom, not just the freely specihable data. 
Second, we could parametrize the gauge invariant content of a general linearized solution in 
terms of 4 functions and on X"*" that capture these true (phase space) degrees of 
freedom. Finally, we showed that the additional condition at sets = 0. Therefore, 
in the linear approximation one sees explicitly that this condition cuts the true degrees of 
freedom in gravitational waves exactly by half. Furthermore, the gravitational waves that 
do satisfy this condition carry no energy-momentum or angular momentum! Thus, although 
this strategy of gaining control over the nonlinear theory seems plausible at hrst, it is simply 
not viable. By isolating the true degrees of freedom at X+, it should be possible to show 
that this sharper results holds also in full general relativity with positive A. 

Second, we found that the conserved energy has a peculiar feature: For matter helds 
as well as linearized gravity, energy ZCt dehned by the time translation can have either 
sign and, furthermore, is unbounded below. Thus, there exist both electromagnetic and 
gravitational waves on de Sitter space-time which carry arbitrarily negative energy, no matter 
how small the positive A is! This is in striking contrast with the A = 0 situation, where 
the corresponding waves carry strictly positive energy Bt in Minkowski space-times. How 
can one reconcile this strong contrast? What happens to the inhnitely many solutions with 
large negative energy in the limit A —)■ 0? 

To analyze this issue, let us hrst recall that, to take this limit, one has to use the differ¬ 
ential structure induced by the coordinates (X, t). In this chart, the cosmological horizons 
which bound region I lie at = (3/A)e“^^* (where = x ■ x). Therefore, in the limit 
A —)■ 0, region I in which is time-like hlls out the whole Minkowski space. This is the 
geometric reason why even though Bt is unbounded below no matter how small the positive 
A is, the limiting T-Lt is strictly positive. In the phase space language, as A changes, the 
covariant phase space Fqq^, on which the Hamiltonian Bt are dehned, itself changes. In the 
limit, the set of solutions hab on which T-Lt is negative simply disappears! 

there are families of metric 


To summarize, as we showed explicitly in section IV B 1 


perturbations 7 a 6 (A) that satisfy the gauge conditions (3.2) and held equations (3.3) for 
each A, and admit well-dehned limits hab as A —)■ 0 satisfying the standard gauge conditions 
and held equations (4.13) in Minkowski space-time. This limiting procedure is onto: the 
limits hab span the entire phase space Fcov of uietric perturbations in Minkowski space. 
Furthermore along any of these families, the energy ZLt 17(A) tends to the energy T-It \f^ of 
the limiting perturbation in Minkowski space. Nonetheless, the lower bound of the energy 
function on phase spaces F^^| is discontinuous in the limit: It equals — cxd for every A, 
however small, but vanishes for A = 0. 

Even though we do recover positivity of energy in the limit A —)• 0, we are left with a 
conundrum because there is strong evidence that A is small but non-zero in our universe: Can 
realistic gravitational waves have arbitrarily large negative energy in de Sitter space-time 
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X+ 


X+ 




FIG. 2: Left Panel: Integral curves of the time translation Killing field T“ in the Poincare patch. 
T“ is future directed and time like in region I and space-like in region II. It is future directed 
and null on portion of the event horizon to the future of the cross-over 2-sphere (bifurcate 

horizon) C and on the portion of the null event horizon to the past of C. It is past directed 

and null on the portion of E~^{i~) to the future of C. Right Panel: A time changing quadrupole 
emitting gravitational waves. Radiation crossing the cosmological horizon E~[i'^) and reaching 
X+ originates in region I; there is no incoming radiation on E~^{i~). 


or, in the nonlinear context, in asymptotically de Sitter space-times? To probe this issue 
let us hrst analyze in some detail the origin of negative energy. Let us begin with Maxwell 
helds in de Sitter space-time. The stress-energy satishes the dominant energy condition and 
the Killing held T“ is future pointing on the part of E^{i~) that lies in region I and past 
pointing on the part that lies in region II (see the left panel in Fig. |^. Therefore, the 
energy hux across E^{i~) into region I is positive but that into region II is negative. It is 
because of this negative hux into region II that the total energy can be negative. Therefore, 
if the Maxwell held under consideration vanished on the part of the horizon E~^{i~) that lies 
in region II, the energy of those electromagnetic waves would be necessarily positive. For 
gravitational waves, we do not have a stress-energy tensor. However, using the fact that the 
Killing held T“ is future directed and time-like in region I, it is easy to show that, if the 
initial data on any cosmological slice S were restricted to he entirely in the intersection of 
S with region I, the energy (4.21) of that cosmological perturbation is necessarily positive.® 
In the limit r] —)■ — oo, the cosmological slice tends to E~^{i~). Therefore, it again follows 
that the conserved energy hux at can be negative only because there is a negative energy 
hux into the Poincare patch across the part of E~^{i~) that lies in region II. But in realistic 
situations gravitational waves from isolated systems would be generated entirely by a time 


® This is most easily seen by using the symplectic form in (4.1) to rewrite the energy as follows: 

1 


Rt = ~ 


SnlLg 


J d^a; ((f"* - 77™) (f” - 77”) 


r , 


+ S™” + 2(1 + -)r^g^ ) 
77 


(5.1) 


where qabf°'f^ = 1 and s™” = g™” — f™f”. If the initial data is restricted to the intersection of E with 
region I we have h/?7| < 1 and, consequently, Ht is necessarily positive. 








20 


changing quadrupole moment (depicted in the right panel of Fig. [^, whence there would be 
no incoming flux across at all. The flux across would just equal that across the 

future horizon that separates regions I and II. Since the Killing held T“ is null and 

future directed on this horizon, this hux has to be positive. Indeed we will show this explicitly 
in m- Thus, in terms of helds at I"^, while general initial data can have arbitrarily large 
negative energies, the initial data induced by gravitational waves produced by realistic sources 
is appropriately constrained for the energy flux across to be positive. In the linearized 
case, it appears to be rather straightforward to make these constraints explicit [15]. An 
interesting challenge in full nonlinear general relativity is to hnd the analogous constraints 
on helds at X+ induced by gravitational waves produced by realistic sources, in absence of 
incoming radiation [T0| (at least from the portion of the event horizon E^{i~) that lies to 
the future of the cross-over surface C). With these constraints at hand, one could hope to 
show that huxes of energy carried by gravitational waves produced in physically realistic 
processes would be positive in full nonlinear general relativity with A > 0, as one physically 
expects. 

Finally, note that our entire analysis -and in particular the limit (4.18) to Minkowski 
space- was carried out by restricting ourselves to the future Poincare patch of de Sitter 
space. As discussed in section [ITA] in the description of isolated systems, this restriction is 
motivated by direct physical considerations. However, one may still ask if the results can 
be extended to full de Sitter space-time. The explicit form of the solutions we presented is 
indeed restricted to the future Poincare patch because of the heavy use of the cosmological 
slicing. But each of these solutions admits a well-dehned extension to full de Sitter space- 
time simply because every solution in our covariant phase space Fcov induces a well-dehned 
initial data on the de Sitter Cauchy surfaces. For these extended solutions, our main results 
also hold on X~. The central formula (4.9) holds for all ten Killing helds of de Sitter 
space-time in this extension. 

These constructions and results provide further guidance for the development of the 
gravitational radiation theory in full nonlinear general relativity with A > 0. We will 
conclude this discussion with two examples. 

Consider hrst the problem of dehning the 2-sphere energy-momentum and angular mo¬ 
mentum charge integrals, analogous to the Bondi 4-momentum in the A = 0 case. For a 
given prescription for selecting asymptotic symmetries, considerations involving held equa¬ 
tions and geometry of X+ (discussed in section 5 of |5]) suggest a natural, candidate expres¬ 
sion for these charges for A > 0. The diherence between these integrals evaluated on any 
two 2-spheres on provides a candidate expression of huxes in the full theory across the 
region of X+ bounded by these 2-spheres. One can show that their linearization provides 
precisely the hux formulas (4.11) and (4.22) at X”^, derived using completely independent 
Hamiltonian methods. This result provides a powerful hint for the charge integrals in the 
full nonlinear theory. The remaining open issue is the selection of appropriate asymptotic 
asymptotic symmetries, without assuming conformal flatness of the intrinsic geometry ofX^ 
(which, as we showed, trivializes the situation by forcing all huxes to vanish). 

A second issue in the full theory is the following. While observations strongly suggest that 
A is positive in our universe, almost all analytical calculations and numerical simulations 
in gravitational wave science set A to zero and work in the asymptotically Minkowskian 
context. (For notable exceptions, see [TGHTS] .) Since the actual value of A is so small 
compared to the scales involved, say, in binary coalescences of astrophysical interest, it is 
natural to assume that setting A to zero is an excellent approximation. However, it is 







21 


not completely clear that this is true for two reasons. First, as we pointed out, the limit 
A —)■ 0 is discontinuous in important respects. Second, advanced LIGO will be eventually 
capable of detecting gravitational waves from sources that are ~ 1 Gpc away, a distance 
that is approximately 20% of the cosmological radius. Therefore, apart from the intrinsic 
conceptual interest, it is important to be able to reliably calculate the ‘errors’ one makes by 
setting A to zero.® The details of the discussion of the A —)■ 0 limit presented in this paper 
will help signihcantly in streamlining these calculations. 
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